Controlled Split-Recombination of 2D Matter- Wave Solitons in Time-Dependent Trap 
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We propose a novel approach to manipulate two-dimensional bright matter- wave solitons by tuning 
the frequency of the trap which is different from Feshbach resonance technique. The exact bright 
soliton solutions for two-dimensional Gross-Pitaevskii (GP) equation with attractive interaction 
strength in a time-dependent trap are constructed analytically and its dynamics show no collapse 
while modulating the trap frequency. The two-soliton dynamics exhibits an interesting splitting 
and recombination phenomenon which generates interference pattern in the process. This type of 
behaviour in two-dimensional BECs has wider ramifications and our approach opens new avenues 
in stabilizing bright solitons in higher dimensional regime. We have also explored the experimental 
realization of this novel phenomenon. 
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I. INTRODUCTION 

A Bose-Einstein condensate (BEC) is a state of mat- 
ter in which a large number of bosonic atoms occupy 
the lowest quantum state of the external potential at ex- 
tremely low temperatures near absolute zero, allowing 
quantum effects to be observed on a macroscopic scale. 
The experimental realization of BECs in weakly interact- 
ing gases has opened the floodgates in the field of atom 
optics and condensed matter physics. The collective exci- 
tation of matter- waves in BECs like matter-wave solitons 
p]-[I| , periodic waves ^ , shock waves ^ , vortices @ and 
necklaces [8| has generated a lot of interest on the explo- 
ration of the dynamics of BECs from both experimental 
and theoretical perspectives. Even though the experi- 
mental realization of Bose-Einstein condensation in dif- 
ferent atomic systems has been observed, the stability of 
BECs imposes restriction on the observation of dynam- 
ical properties of such topological excitations. This is 
because of the higher dimensionality of the system. For 
example, attractive Li BECs have been shown to collapse 
in three dimensions [H, [loj- But, the one- dimensional at- 
tractive Li BECs in a standing light wave potential are 
stable This is primarily because of the fact that 

in higher dimensional BECs, a nonlinear excitation of 
atoms just bigger than the background is susceptible to 
even smaller perturbations. The dark solitons will decay 
into vortex pairs under small transverse perturbation [7j . 
Therefore, longtime dynamical behaviour of nonlinear ex- 
citations in higher dimensional BECs is hard to observe 
in real experiments. 

However, it was identified that the bright solitons and 
vortex solitons can be stabilized without trap in higher 
dimensional BECs by rapidly oscillating atom-atom in- 
teraction strength using Feshbach resonance technique 
[l2l-[l4j where resonantly tuning uniform magnetic field 
produces an effective confinement creating stable self- 
confined condensates. Similarly, a stable shallow ring 
dark soliton has also been observed in two-dimensional 
BECs with tunable interaction [l5|. Thus, to stabilize 



the solitons in BECs, one has to change the interaction 
strength between the atoms by tuning magnetic field near 
Feshbach resonance [T^, 03 • 

Can one suitably modulate the trap frequency in a 
time-dependent trap and control the stability of the 
condensates instead of varying the interaction strength? 
This paper is aimed at the investigation of the (2+1) 
dimensional GP equation with time-dependent trapping 
potential describing the dynamics of matter waves in 
pancake-shaped BECs. The fact that the trap frequency 
varies with time means that one can suitably tune the 
time-dependent harmonic trap to stabilize the bright soli- 
tons and this is an alternative way of stabilizing the 
bright solitons in BECs rather than tuning the interac- 
tion strength. We first construct exact bright matter- 
wave soliton solutions using Hirota method [18j and then 
study the impact of modulating the trap frequency. We 
demonstrate split-recombination phenomenon by mod- 
ulating the trap frequency. The frequency modulation 
technique is well known in controlling molecular BEC 
scattering length and condensate spilling [19]. Our re- 
sults open a new window to manipulate atoms in higher 
dimensional BECs. 



II. MODEL 

At ultra low temperatures, the dynamics of BECs can 
be described by the three dimensional time-dependent 
GP equation of the following form [20] 



ih 



dt 



2m 



V(v,t) + U\*(v,t)\ 2 ^(v,t) 

(1) 

where ty[r(x, t), t] represents the condensate wave func- 
tion normalized by the particle number N = J <ir|\I/| 2 , V 2 
denotes the Laplacian operator and V(r) is the trapping 
potential of the form V(r) = m(uu r (t) 2 r 2 + uj 2 z 2 ), where 
r 2 = x 2 + y 2 , uOr :Z are the confinement frequencies in the 
radial and axial directions, respectively. The coefficient 
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U is related to s-wave scattering length a s controlled by 
Feshbach resonance as U = 47rft 2 a s /m, m is the atom 
mass. 

The three-dimensional time-dependent GP equation 
effectively reduces to two-dimensional GP equation when 
the radial trapping frequency u r is very small compared 
to axial trapping frequency uj z , i.e (\uj r \/uj z <C 1). Con- 
sequently, the atoms move freely in the x, y plane. Al- 
lowing the radial trapping frequency to vary with time 
i.e. uj r = u) r (t), the corresponding two-dimensional GP 
equation becomes 

1 dt + 2 \dx 2 + dy 2 

where the strength of trapping potential is related to ra- 
dial and axial trapping frequency as Q(t) 2 = \uj r (t)\ 2 / uo 2 
and g = (SirniuUz/hY^Nas represents the strength of 
attractive interaction (g > 0). In Eq. (2) the length, 
time, frequency and wave function are measured in units 
of an, 1/Wo, ^ro, and 7V 1 / 2 /a/ l , respectively, where 
ah = y/h/muj r o and u r o is a radial trapping frequency at 
t — to. When Q(t) 2 > 0, the trapping potential is expul- 
sive and Q(t) 2 < represents the confining potential. 

The nature of the nonlinear excitations and stability 
of pan cake shaped BECs in unmodulated traps (time 
independent traps) is well studied in different physical 
contexts [2~0-[23l]. However, the impact of modulating 
the trap frequency on the condensates is not yet known. 
So, the investigation of (2+1) GP equation in an arbi- 
trary time dependent trapping potential with constant 
attractive interaction strength assumes tremendous sig- 
nificance. 



It should be mentioned that the choice of trapping poten- 
tial strength depends on the solvability of Ricatti equa- 
tion (Eq. (6)), but the solution of Ricatti equation is 
not unique. One could also get different physically real- 
izable potential by properly choosing the time dependent 
function a(t). 

The operator D x and D t are defined as 



D?D™{G ■ F) 



d d 



d d 



dt dt' J \dx dx' 



G(t, x)F(t', x')\ t '=t,x'=x- 



ip + g\ijj\ 2 ip + = 0, (2) The functions G and F can be extended as 

G = £5 (1) +e¥ 3) +£ 5 5 (5) + -- - , 

F= l +e 2 / (2) +£ 4 / (4) + ... _ 



(7) 



(8) 
(9) 



Substituting G and F into the bilinear forms and col- 
lecting the coefficient of various powers of e, one obtains 
the system of linear partial differential equations which 
can be recursively solved. The function g^ can be rep- 
resented as the plane wave solution of the form 

N 

gd) = J2 e * , Xj = h[ j) (t)x + (t)y + hf (t). 

3 = 1 

(10) 

To generate one soliton solution, we take N = 1 to 
obtain g^ and as 



^(1) =e Xi(x,y,t)^ 

y(2) = e xi(x tyi t)+xU^y^)+m(t) 



(11) 
(12) 



III. ONE SOLITON SOLUTION 

The aim of this paper is to solve the (2+1) GP equation 
for an arbitrary time-dependent trapping potential with 
constant attractive interaction strength. To solve Eq. 
(2), we use Hirota bilinear method. To generate Hirota 
bilinear form of Eq. (2), we make use of the following 
transformation of the form 



with 



i/j = e 



F' 



(3) 



where G = G(x,y,i) is complex, while F = F(x,y,t) is 
real. Substituting Eq. (3) into Eq. (2), one obtains the 
following Hirota bilinear form, 

(iD t + ^{D 2 x + D 2 y ) + ia(t)xD x +ia(t)yD y 

+ ia(t))G • F = 0, (4) 
\{Dl + Dl)F.F = g\G\\ (5) 

and the trapping potential f£(t) 2 is related to the function 
a(t) as 



90^2 da ( f ) | „v^2 



(6) 



Xi(x,y,t) = h^(t)x + h^(t)y + h^(t), (13) 
hf\t) =ae-f a ^ dt , h^(t) = be-f a ^ dt , 

h ( s\t) = JC-(a 2 +b 2 )e- 2 f^ dt -a(t))dt,(U) 

e 2fa(t)dt 



(a + a*) 2 + (b + b*) 2, 



(15) 



where a and b are complex constants, and xt denotes 
complex conjugate of xi? is real function. Substi- 
tuting g^ and into the truncated series, we obtain 
G = and F = 1 + f^ 2 \ Thus, we get the one soliton 
solution of the following form 



l -e-^ 2 sech (iJe(xi) + y 



JIm( X i)+ia(t)/2(x 2 +y 2 ) _ 

(16) 

From the above solution, we observe that the ampli- 
tude (e- 2 f a ^ dt (a 2 R + b^/Ag) 1 ' 2 , velocity (Re (h 3 (t)) + 
rji(t)/2), position(Re(xi) + ?7i/2) and phase (arg(^)) of 
the soliton depend on the frequency of the trapping po- 
tential (i.e. 2u r (t) 2 /u 2 z = 2Vt(t) 2 = a'{t) + a(t) 2 ) and 
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attractive interaction strength g. Since the frequency of 
the trapping potential could be modulated with respect 
to time, one can easily control the physical properties of 
bright solitons and hence the condensates. We can choose 
the trapping potential to be either confining (Q(t) 2 < 0) 
or expulsive (Q(t) 2 > 0). And, we can also continu- 
ously change the nature of the trapping potential from 
confining to expulsive or expulsive to confining. It is in- 
teresting to note that this exact bright soliton solution 
involves constant attractive interaction strength which 
is independent of trapping potential strength unlike the 
quasi one-dimensional BECs [2J-[26|, wherein one needs 
a delicate balance between time-dependent trap and in- 
teraction strength. 

The time-dependent trapping frequency allows us to 
choose different physically realizable potentials and study 
the effect of trap on the dynamics of solitons. For 
example, in the case of a constant trapping potential 
(Q(t) 2 = Qo — 0.4), the amplitude and width of the 
soliton varies with time as shown in Fig. 1(a), in which 
the amplitude of the soliton increases and reach the max- 
imum value which depends on the strength of trapping 
potential and then decreases. When the trapping po- 
tential is suddenly switched off, the amplitude decreases 
and the width increases as shown in Fig. 1(b). From 
these two cases, we observe that the time-independent 
trap and zero-trapping potential cannot keep the soli- 
ton as stable. As an alternative, one can look to stabilize 
the bright solitons by modulating the trapping frequency 
with time. For example, choosing Q(t) 2 = — 2 — cos(2t), 
the a(t) takes the following form 



a(t) 



mathieuC '(4, —1, t) + mathieuS' '(4, —1, t) 
mathieuC(A, —1, t) + mathieuS(A, —1, t) 



The mathieuC (Z, m, t) and mathieuS (/, m, t) are even and 
odd functions, respectively. When m = 0, they are sim- 
ply as cos(Vlt) and sin(V^) where the prime functions 
are t derivatives of matheiu functions. When we modu- 
late the trap strength periodically by tuning the radial 
trapping frequency (uj r (t)) within the two-dimensional 
confining regime, the amplitude and width of the soliton 
oscillates and gets amplified as shown in Fig. 1(c), but 
does not collapse or dilute even if one waits long enough. 
In this case, the stability of the soliton is partial. One 
can choose the other possible time-dependent potentials 
to stabilize the soliton where the amplitude and width of 
the soliton remains constant or oscillates without ampli- 
fication. 

When the confining trapping potential strength de- 
creases exponentially (ft(£) 2 = e~ 4£ (0.1352 - 0.52e 2£ )) 
as shown in Fig. 1(d), the amplitude and width be- 
comes constant for large t. This means that the sta- 
bility of the soliton is ensured. If we choose the trap- 
ping strength to oscillate between confining and expulsive 
regime (Vt(t) 2 = -0.5cos(2£) + 0.125 sin(2t) 2 ) as shown 
in Fig. 1(e), both the amplitude and width of the soli- 
ton oscillate but no amplification occurs. When we in- 
crease the frequency of oscillations of trapping frequency 
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FIG. 1: (Color online) The amplitude, width and trap- 
ping strength of one-soliton solution plotted for different 
case, (a) For constant time-independent trapping potential 
(Q(t) 2 = 0.4 (i.e., Ur = 2?r x 284Hz, u z = 2tt x 710Hz)) with 
g = 0.0005, a = 10 + 0.2z, b = 10 - 0.3z. (b) The trapping 
potential is switched off when g — 0.05, a = 1.5 + 0.2i,b = 
1.5 — 0.3i. (c) Oscillating time-dependent trap within con- 
fining regime (^(t) 2 = —2 — cos(2t)), (d) exponentially de- 
creasing confining trap approaching zero, (e) trapping fre- 
quency oscillates between confining and expulsive regime and 
(f) increasing oscillation frequency of the trap with g = 0.05, 
a = 0.15 + 0.2i,6 = 0.15-0.3i. 



(tt(t) 2 = -2cos(8t) + 0.125 sin(8t) 2 ) as shown in Fig. 
1(f), the amplitude and width oscillates but the oscil- 
lation is very small. In other words, the frequency of 
oscillation of the amplitude is proportional to frequency 
of oscillation of time-dependent frequency. This is sim- 
ilar to the Feshbach resonance [12| technique in which 
the amplitude of oscillating soliton is nearly a constant 
while increasing the frequency of interaction strength. 
From the Figs, l(a-f), we observe that the longevity of 
soliton can be controlled by fine tuning the trapping po- 
tential strength. One can realize the above phenomenon 
in experiments as well. For example, in the case of 7 Li, 
keeping BEC the attractive interaction strength g = 0.25 
and uj z = 2tt x 710Hz, one can oscillate the radial trap- 
ping frequency approximately between 2tt x 80Hz and 
2tt x 200Hz. Invoking this ideas in the two-soliton solu- 
tion, one observes an interesting splitting and recombi- 
nation phenomenon. 
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IV. SPLIT AND RECOMBINCATION O 
TWO-SOLITON 

To generate two-soliton solution, we consider N - 



the expression = ^2f =1 e Xj in Eq. (10) and r 
the same procedure as in one soliton case, to get 
9^\ an d while the series gets truncated al 
The two-soliton solution can be explicitly w 

as 



^2 
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(i) 
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(3) 



1 + /(2) + /(4) ; 



where 



^(3) 
K2) 



K4) 



e xi+xi+X2+ri _^ e xi+X2+x2+ r 2 ^ 

e Xi+xI+^7i _|_ e X2+xS+^72 
gXi+X^+^o ^ e xt+X2+^7o^ 



and 



Xi 


= hf (t)x + hf {t)y + hf (t) 


; i = i,: 






= ae-f a ^ dt ; h% ] = be~ ^ 


i(t)dt 
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h ( s\t) 
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a(t))dt, 




e r > 


= fte" 1 +<5 2 e"S, e r2 = <5 3 e" 


+ (5 4 e" 2 , 




e" 1 


= < ? e 2 J ra(t)dt /((a + o*) 2 + (6 4 


-6*) 2 ), 




e" 2 


= 5 e 2 i' a W d V((c + c*) 2 + (rf + 


-d*) 2 ), 


(29) 




= f/e 2 J ra( * )d 7((a + c*) 2 + (6 4 


-d*) 2 ), 


(30) 




= 5 e 2 i' Q W d V((c + a*) 2 + (c/4 


-b*)\ 


(31) 



h = 

S 3 = 

6 4 = 

m(t) = 



(c~ 


- a)(a H 


-a*) 


+ (d-b)(b + 


b*) 


(c + 


a*) (a - 


-a*) 


4(d+6*)(6- 


f 6*)' 


(a 


-c)(c4 


-a*) 


+-(6-d)(d + 


6*) 


(c + 


a*) (a 4 


-a*) 


+ (d + 6*)(6- 


f &*)' 


(c- 


- a) (a - 


-c*) 


4 (d -&)(& + 


d*) 


(c + 


c*)(a 4 


-c*) - 


f (d + d*)(&- 


fd*)' 


(a 


-c)(c4 


-c*)- 


f (6-d)(d + 


<f) 


(c + 


c*)(a4 


-c*) - 


f (d+d*)(b- 


f d*)' 


1 




(t)dt( 


a — c) 2 (a* - 


c*) 2 


4 (a 


4a*) 2 ( 


c + a 


*) 2 (a4c*) 2 ( 


c4c*) 2 ) 



(32) 
(33) 
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(35) 

(36) 



The constraint a = b and c = d ensures the two-soliton 
solution to be exact and the solitons are parallel to each 
other so called "parallel bright matter- wave solitons" . 
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FIG. 2: (Color online) (a) The split-recombination of two par- 
allel bright matter-wave solitons with the following parame- 
ters a = b = 0.3 + 0.H, c = d = 0.3, a(t) -0.001(t - 50). 
(b) The snapshots of two-soliton along z-axis at different in- 
stants of time. The trap strength Q(t) 2 parabolically varying 
between —0.0001 and —0.0005 in a span of 80 ms. 



From the two-soliton solution (Eq. (17)), one can ob- 
serve the dynamics of parallel bright matter-wave soli- 
tons in the two-dimensional harmonic trap. We choose 
the function a(t) = — 0.001 (t — 50) in which the trap 
strength (2Vt(t) 2 = a'(t) + a(t) 2 ) varies within the con- 
fining regime (Q(t) 2 < 0) between t ~ 10 and t ~ 80. 
First, we consider a two parallel bright solitons merged 
in a single bound state with the initial conditions a = 
b = 0.3 + 0.H,c = d = 0.3, and attractive interaction 
strength is kept constant at g = 0.5. We then modu- 
late the trap strength in a parabolic manner as shown 
in Fig. 2(a). From Fig. 2(a), we observe that a soli- 
ton coherent state at t = 10 evolves in time and splits 
into two parallel solitons. Later on, they recombine at 
t = 80 and produce an interference pattern as shown in 
Fig. 2(b). The formation of interference pattern due to 
the collision of two parallel bright solitons with zero angle 
of incidence in a single component condensate is different 
from the interference pattern produced by the expansion 
of two condensate clouds [27j when the trapping poten- 
tial is switched off. Thus, our results reinforce the fact 
that the matter-waves originating from the condensates 
(bright solitons) do interfere and produce a fringe pattern 
as a clear signature of the long range spatial coherence of 
the condensates. This is another way of generating inter- 
ference pattern and one can delay or control the creation 
of pattern by properly modulating the trapping potential 
strength. The two parallel bright soliton maximum sep- 



5 



aration depends on the strength of trapping potential as 
shown in Fig. 2(a). Experimentally, one can realize this 
phenomenon by parabolically modulating the trapping 
frequency within the two-dimensional confining regime. 

In experiment, to generate parallel bright solitons in 
a two-dimensional harmonic trapping potential, initially 
the trapping frequencies in the case of 7 Li BECs are 
uj r = 2tt x 50Hz and uo z = 2tt x 710Hz with the effective at- 
tractive interaction being g = 0.25. This trap can be de- 
termined by a combination of spectroscopic observations, 
direct magnetic field measurement and the observed spa- 
tial cylindrical symmetry of the trapped atom cloud [28| . 
After making this set up to generate two-soliton, the con- 
densates split into two parts using radio frequency (Re- 
dressed potential [29]. This technique has the advantage 
of allowing a smooth transition from a single trap into 
a double well potential and hence coherent splitting is 
possible with splitting range from 3 — 80/im. Due to 
the constant attractive scattering length, the two-bright 
solitons are generated from two nearly separated conden- 
sates. We may consider this as a initial coherent state of 
two-soliton split-recombination phenomenon. Then, one 
starts to tune the radial trap frequency in a parabolic 
manner after switching off the RF. When the radial trap- 
ping frequency which after fine tuning reaches the ini- 
tial value, the two-bright solitons recombine and produce 
fringe pattern. 



V. CONCLUSION 

We have found exact bright matter-wave solitons 
for a two-dimensional BEC with attractive interaction 



strength in a time-dependent harmonic trap. Even 
though attractive interaction strengths usually lead to 
the collapse of the condensates, one can sustain the sta- 
bility of solitons by properly modulating the frequency 
of the radial trap. Our method is different from stabi- 
lizing solitons by Feshbach resonance technique. We also 
observe that the splitting and recombination of matter- 
wave solitons which generates interference pattern in the 
process is a new phenomenon in two dimensional BECs 
and may have wider ramifications in the manipulation 
of atoms in BECs. In particular, ultra-narrow two- 
dimensional BEC solitons is very useful in the field of 
nanolithography. 
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